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PREFACE 


In  a  previous  report  (AMP  Report  42*11),  the  problem 
of  vertical  Impact  of  a  body  with  a  spherical  nose  on  a 
water-surface  was  solved  mathematically  by  Max  Shiffman  arx3 
Donald  C*  Spencer  of  the  Few  York  University  Group  of  the 
Applied  Mathematics  Panel*  The  solution  was  based  on  the 
assumption  that  the  surface  of  the  water  is  not  appreciably 
altered  by  the  entry  of  the  body,  an  assumption  which  can 
be  accepted  only  as  a  first  approximation* 

Therefore,  a  further  study  was  undertaken,  taking 
into  account  the  disturbance  of  the  surface  of  the  water* 

The  present  report  presents  the  results  of  these  studies* 

It  is  shown  that  the  disturbance  of  the  surface  leads  to  two 
types  of  corrections  of  the  previous  results*  In  the  first 
place,  the  upward  rise  of  the  surface  of  the  water  wets  a 
larger  portion  of  the  entering  body,  thereby  increasing  the 
resistance  of  the  water  to  the  entering  body*  Secondly,  the 
rising  surface  relaxes  somewhat  the  restraint  imposed  on  the 
body,  causing  a  decrease  in  the  impact  force.  These  two 
corrections,  which  partly  counterbalance  one  another,  are 
analyzed  here  by  using  a  combination  of  theory  and  information 
derived  from  experiment*  The  final  result,  which  deviates 
somewhat  from  the  first  approximation.  Is  compared  with 
experiment  and  shown  to  be  in  satisfactory  agreement  with  it* 

With  the  experience  gained  and  the  theory  developed  for 
vertical  entry,  it  is  hoped  that  the  impact  forces  during 
oblique  entry  will  be  amenable  to  treatment*  This  question 
Is  now  being  considered  by  the  same  authors* 

Attention  should  be  drawn  to  the  work  of  the  Harvard 
Group  of  the  Applied  Mathematics  Panel  concerning  air-water 
entry  as  a  whole,  and  to  the  activities  of  the  Underwater 
Ballistics  Committee* 

The  writers  were  greatly  helped  by  other  members  of  the 
Hew  York  University  Group,  in  particular  by  Mr*  M.  H.  Shamos, 
who  niade  an  important  contribution  by  excellent  flash  photo¬ 
graphs  he  took  of  spheres  entering  water.  These  photographs 
furnished  Bruch  qualitative  and  numerical  information*  Thanks 
are  also  due  to  Dr.  Patrick  Hurley  of  the  Morris  Dam  Group  of 
the  California  Institute  of  Technology  for  permission  to  use 
the  photographs  exhibited  Id  Plate  5;  to  Dr.  R.  M.  Davies  of 
the  Engineering  Laboratory,  Cambridge,  England  for  Plate  6; 
and  to  Dr.  E.  Newton  Harvey  of  Princeton  University  for  Plate  7. 

R.  Courant 

Director,  Contract  OEMsr-945 
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THE  FORCE  OF  IMPACT  ON  A  3PHERE  STRIKING  A  WATER  SURFACE 

Second  Approximation 

1.  Introduction  and  Summary 


The  entry  of  a  solid  "body  from  air  into  water  may 
take  place  in  two  ways:  (i)  if  the  speed  is  small  and  if 
the  surface  of  the  body  is  smooth,  the  entry  occurs  with¬ 
out  the  formation  of  a  cavity  (smooth  entry);  (11)  if  the 
speed  is  great  or  if  the  surface  of  the  body  is  rough  or 
has  an  irregular  shape,  the  flow  detaches  from  the  body 
when  it  has  penetrated  a  short  distance  below  the  initial 
surface  level  (rough  entry).  In  either  case  there  is  an 
initial  period  following  contact  with  the  surface  during 
which  the  body  experiences  the  greatest  deceleration  and 
which  is  of  such  short  duration  that  cavitation  has  not 
yet  developed.-^  The  determination  of  the  forces  acting 
during  this  impact  stage  of  the  entry  is  important  in  the 
case  of  pro^otiles  shot  or  dropped  into  water,  not  only 
because  oar  possible  damage  to  the  mechanism  and  nose  but 
also  because  the  impulsive  forces  and  torques  created' by 
the/impact  influence  the  underwater  trajectory  of  the  pro- 
>ctile.  1  , ,  ,w  :'/  L 

Wo  cu.ppc.se^ that  the  nose  of  the  projectile  is  spheri¬ 
cal  in  shape  and.that  the  initial  inpact  with  the  surface  is 
vertical. 


then  visualizeithe  mass  of  the  projectile'* 
as  concentrated  into  a  sphere  of  the  same  radius  as  the  nose. 
In  a  previous  r ep or t\ "ob t a lnedja  first  approximation; 
to  the  foroe  of  impact  on  a  sphere  entering  vertically  by 
disregarding  the  rise  of  the  surface;  in  the  present  report 
a  second  approximation^ by "estimating  the  effects 

\Of  the  surface  motion  and  .i 

X  O  -vA 


a.  with  experiments^.- 


-ouar  theoretical  estimates  /U\ 


*  Humber s  in  bra 
of  this  report 


_ts  [  ]  refep  o  the  bibliography  at  the  end 
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The  resultant  upward  force  P  acting  on  the  sphere 
In  the  Impact  stage  can  he  expressed  In  terms  of  a  dimen- 
slonless  Impact-drag  coefficient  C  •  Let  A  be  the  radius 
of  the  sphere,  TJ0  its  Initial  velocity  upon  Impact,  and 
p  the  density  of  water#  Then  the  impact  force  may  be 
conveniently  expressed  as 

p  cp 

where  Cp  is  the  "impact -drag  coefficient# "  The  coefficient 
C  depends  on  the  time;  it  is  zero  when  the  sphere  first 
strikes  the  surface,  rapidly  increases  to  a  maximum  value, 
and  then  decreases# 


The  impact-drag  coefficient  Cp  may  be  expressed  in 
terms  of  the  "virtual  mass"  of  the  fluid#  If  XT  is  the 
velooity  of  the  sphere  at  time  t,  we  define  the  virtual 
mass  of  the  fluid  to  be  a  quantity  M  such  that  MU  Is  the 
total  vertical  inpulse  communicated  to  the  fluid  during 
the  time  interval  t;_  in  other  words  M  is  defined  by  the 
equation 

p  =  <Mt7)  • 

Let  B  denote  the  depth  of  penetration  of  the  bottom  of  the 
sphere  below  the  initial  surface  level,  and  let  b  =  B/A# 

mmsmsaam 
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1  3 

Introducing  the  additional  dimensionless  quantity  m  =s  M/g  VA* 
we  obtain  In  Section  3  the  .formula 


/ 


°p  = 


dm 


where  <r  is  the  "specific  gravity"  of  the  incoming  projectile 
considered  as  though  the  whole  mass  were  concentrated  into 
a  sphere  having  the  sane  radius  as  the  nose, 

_  _  mass  of  projectile 

(j  a  1  11  "  ■■■"  —  -  • 

fv*5 

Our  final  estimate  of  Cp  as  a  function  of  b  is  plotted 
in  Graph  1  for  various  valuos  of  «■  .  (All  graphs  are  collected 
together  at  the  end  of  the  report*  For  the  sake  of  convenience. 
Graph  1  is  also  inserted  on  the  next  page*)  For  values  of  <r 
ranging  from  1  to  oo  ,  the  impact  drag  coefficient  Cp  rises  very 
rapidly  to  a  peak  of  about  1,  the  peak  occurring  when  the  sphere 
has  penetrated  to  a  distance  of  from  0*1  to  0*2  of  a  radius 
below  the  initial  surface  level*  The  coefficient  Cp  then  de¬ 
clines  more  slowly,  until  at  about  0.7  of  a  radius  of  entry 
its  value  is  equal  to  the  stationary  cavity  drag  value  of  0*3 
(far  description  of  cavity  drag,  see  [2]  and  [21]). 

For  smaller  values  of  cr  ,  this  general  behavior  of  Cp 
is  retained  but  is  shifted  downward  and  to  the  left* 

The  question  arises  whether,  in  rough  entry,  the  cavity 
formation  which  takes  place  will  interfere  with  the  above  re¬ 
sults  and  partly  vitiate  them.  It  is  known  from  both  experimen¬ 
tal  and  theoretical  grounds  that  the  cavity  begins  to  form^fely 
after  the  sphere  has  submerged  a  distance  somewhere  between  a 
quarter  of  a  radius  and  a  radius.  Thus  the  behavior  of  Cp  de¬ 
picted  in  Graph  1  still  applies  —  the  rapid  rise  of  Cp  to  a 
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maximum  of  about  1,  and  the  slower  decline  to  the  cavity 
drag  value  of  0.3.  The  brief  impact  stage  in  rough  entry 
is  the  same  as  in  smooth  entry. 

The  dimensionless  virtual  mass  m  defined  by  m  =  VL/±is  pk 
is  plotted  as  a  function  of  b  in  Graph  3*  For  purposes  of  com¬ 
parison  we  have  plotted  in  the  same  graph  our  first  approxima¬ 
tion  for  m,  taken  from  the  previous  report  [l5]«  Likewise, 
in  Graph  2,  our  corrected  estimate  of  Cp  in  the  case  <r  =  oo 
is  compared  with  our  first  approximation  obtained  in  [15].  In 
the  corrected  estimate,  the  maximum  of  Cp  is  somewhat  larger 
and  Is  displaced  to  the  left. 

The  first  approximation  of  the  Impact-drag  coefficient 
was  derived  in  our  report  C 15]  under  the  assumption  that  the 
surface  does  not  rise,  the  surface  remaining  plane  and  at  its 
initial  level.  Under  these  circumstances  the  flow  at  each 
instant  t  depends  only  on  the  depth  of  penetration  below  the 
initial  level  and  on  the  velocity  of  the  sphere  at  that  Instant 
t  and  Is  otherwise  independent  of  the  past  history  of  the  entry; 
for  this  reason  the  flow  might  be  called  " quasi -steady. n 

In  the  present  report  a  second  approximation  Is  obtained 
which  takes  into  account  the  past  history  of  the  phenomenon. 

The  corrections  are  of  two  sorts.  On  the  one  hand  the  water 
rises  and  wets  a  larger  portion  of  the  spherical  surface  than 
it  would  if  the  surface  did  not  rise.  The  fluid  therefore 
exerts  pressure  over  a  larger  area  and  Increases  the  total 
impact  force.  This  will  be  called  the  "wetting  correction.* 

On  the  other  hand,  the  free  surface  rises  and  relaxes  to  a 
certain  extent  the  restraint  imposed  on  the  sphere.  This  will 
be  called  the  "surface  correction."  These  two  corrections 
are  discussed  in  Sections  6  and  7  respectively,  and  both  are 
combined  in  Section  8  to  give  the  final  estimate  of  the  inpact- 
drag  coefficient  Cp. 

In  Section  9  this  final  estimate  of  Cp  is  compared  with 
experimental  measurements  of  the  Impact  force  on  a  sphere. 


Comparison  is  mads  with  measurements  taken  'oy  S.  Watanabe’s 
[17]  In  1930,  who  used  a  piezoeleotrlo  gauge  to  determine 
the  force;  with  rough  experiments  of  R.  W.  Blundell  [3]  made 
in  1937;  and  with  recent  measurements  of  E.  G.  Richardson  [  13]  • 
Watanabe’s  experiments  were  very  carefully  performed 
with  accurate  instruments  and  verified  the  fact  that  the  impact 
force  varies  as  the  square  of  the  initial  velocity.  The  value 
of  the  spocifio  gravity  <r  of  the  sphere  in  Watanabe's  experi¬ 
ment  was  as  low  as  0*116.  The  theoretical  Cp  curve  correspond¬ 
ing  to  this  value  of  <r  is  compared  with  Watanabe’s  experimental 
curve  in  Graph  10*  Vote  that  the  range  of  b  is  small,  from  0 
to  *03.  The  agreement  is  remarkably  good  for  very  small  values 
of  b,  when  Cp  is  increasing  rapidly.  In  fact,  near  b  =  0  the 
theory  gives 

Cp  »  6.62  b1^2  -  a2b  «  ... 

where  the  constant  is  not  accurately  determined  by  the  theory. 
The  first  term  6.62b  of  this  expansion  is  verified  completely 
by  Watanabe’s  data*  The  deviation  in  the  remaining  portion  of 
the  curve  is  due  to  the  second  coefficient  a2* 

In  Blundell »o  experiment,  <r  ranged  from  .2  and  .3  approxi 
mately.  The  experimental  curves,  transcribed  to  the  dimension¬ 
less  form  of  a  drag  coefficient,  are  shown  in  Graph  11(a),  and, 
the  corresponding  theoretical  curves  are  drawn  in  Graph  11(b). 

An  idea  of  the  rather  large  experimental  error  may  be  seen  by 
comparing  curves  (B)  and  (C)  which  should  be  identical,  and  from 
the  fact  that  curve  (D)  should  be  the  highest  (since  <r  is  the 
largest).  The  experiments  at  least  give  an  idea  of  the  order  of 
magnitude  of  the  impact  ibrces  and  of  their  duration,  and  the 
agreement  with  theory  Is  satisfactory.  In  fact,  curve  (B)  even 
agrees  numerically  with  theory. 

In  Richardson’s  experiment,  <r  -  .16.  The  corresponding 
theoretical  curve  Is  compared  with  Richardson’ s  in  Graph  12, 


and  the  agreement  is  fair*  Again,  Richardson*s  experiment 
should  he  considered  as  yielding  only  the  order  of  magnitude 
of  0p. 

In  each  of  these  experiments,  the  value  of  <T  is  small* 

It  would  be  desirable  to  carry  out  experiments  with  larger 
values  of  <r  ,  using  accurate  instruments  (such  as  piezoelectric 
gauges)* 

In  addition  to  the  inpact  force,  it  is  of  interest  to 
know  the  average  pressure  p  acting  over  the  portion  of  the 
sphere  in  contact  with  the  water*  This  pressure  can  be  written 
in  the  form 

P  =  *  Dp 

where  Dp  is  a  dimensionless  coefficient*  The  dependence  of 
Dp  on  b  =  B/A  is  shown  in  Graph  4  (for  the  case  <r  =  oo  ),  it 
is  noteworthy  that  the  value  of  Dp  at  tin  initial  contact  of 
the  sphere  with  the  water  is  infinite*  Actually,  because  of 
the  compressibility  of  the  water  and  because  of  the  fact  that 
the  compression  wave  is  initially  a  plane  wave,  the  pressure  p 
during  this  phase  is  approximately  cU0  where  c  is  the  speed 
of  sound  in  water*  This  would  give  a  value  of  Dp  approximately 
equal  to  •  The  intersection  of  the  straight  line  Dp  =  ^ 
with  Graph  4  gives  an  indication  of  the  time  interval  in  which 
compressibility  plays  a  role*  Thus,  the  effect  of  compressibility 

*  Using  the  asymptotic  formula 


for  small  b,  and  equating  this  to  -£■  ,  one  obtains 

o  U° 

as  an  estimate  of  the  relative  depth  of  penetration  in  which 
compressibility  plays  a  role*  Even  if  the  initial  entry 
velocity  UQ  is  as  high  as  500  ft/seo  this  yields  b  =  vOl,  an 
order  of  magnitude  which  is  negligible* 

This  idea  was  also  suggested  by  M*  Gimprich. 
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Is  to  reduce  the  value  of  Dp  and  therefore  of  Cp  at  the 
very  beginning. 

Concerning  the  actual  distribution  of  pressure  over 
the  wetted  portion  of  the  sphere,  we  remark  that  the  pressure 
is  not  greatest  at  the  lowest  point  of  the  sphere  (stagnation 
point).  This  fact  arises  from  the  non -stationary  character 
of  the  flow,  and  the  usual  conceptions  based  on  stationary 
flow  do  not  apply.  On  the  contrary,  the  pressure  is  a  mini¬ 
mum  at  the  lowest  point  of  the  sphere  and  has  a  very  sharp 
and  high  roaxlmnm  near  the  edge  of  the  wetted  portion  of  the 
sphere*  This  fact  was  derived  theoretically  by  Wagner  [16] 
in  the  two-dimensional  case  of  the  entry  of  a  very  flat  wedge. 
However,  an  estimate  of  the  depth  at  which  stationary  flow 
approximately  sets  in  can  be  obtained  by  determining  when  Cp 
Is  about  0*3.  This  oocurs  for  b  *=  *7  approximately. 

Sootlon  2.  Photographs  of  Entry 

The  interesting  sequence  of  events  produced  by  entry 
of  a  sphere  from  air  into  water  was  first  photographed  and 
described  by  Worthington  [19].  Plate  1  is  a  reproduction 
of  some  of  his  photographs  and  shows  several  stages  in  the 
smooth  entry  of  a  sphere*  A  more  couple te  sequence,  which 
show-  the  splash  in  greater  detail,  is  given  in  Plates  2,3. 

As  may  be  observed  from  these  two  sets  of  pictures  the  phenom¬ 
ena  accompanying  smooth  entry  are  closely  duplicated  in  differ¬ 
ent  experiments*  Even  the  Jets  produced  as  the  water  envelopes 
the  spheres  differ  little  in  appearance* 

In  order  to  show  more  clearly  the  splash  during  impact, 
an  enlarged  photograph  is  reproduced  In  Plate  4.  We  observe 
that  the  splash  Is  largely  concentrated  in  the  immediate 
neighborhood  of  the  sphere  and  is  composed  of  a  relatively 
thick  base  from  which  a  thin  sheath  rises  up  along  the  sides 
of  the  sphere.  A  few  drops  which  have  detached  themselves  may 
be  seen  at  the  top  of  the  sheath. 


Because  of  the  speed  with  which  projectiles  are 
usually  launched  into  water  their  entry  is  rough.  In  fact, 
the  entry  of  a  boni),  mine,  or  torpedo  takes  place  in  three 
stages:  (i)  the  Initial  impact  stage  which  begins  when  the 

projectile  touches  the  surface  and  ends  when  cavitation 
developes;  (11)  motion  through  the  water  with  a  trailing 
cavity  extending  to  the  surface;  (iii)  motion  after  the 
cavity  is  closed.  The  closure  of  the  cavity  takes  place 
either  at  the  surface  or  below  the  surface,  but  in  either 
case  the  projectile  carries  with  it  a  pulsating  air  bubble 
which  gradually  disappears. 

The  initial  stage  of  the  rough  entry  of  a  sphere  is 
shown  in  Plate  5.  The  sheath,  which  persists  alongside  the 
sphere  in  the  case  of  smooth  entry,  detaches  in  rough  entry 
and  may  break  up  into  spray.  However,  the  thick  base  of  the 
splash  remains  in  contact  with  the  sphere  and  is  the  same  as 
in  smooth  entry.  So  far  as  estimation  of  impact  force  is 
concerned,  the  effect  of  the  sheath  is  negligible  compared  to 
that  of  the  thick  base. 

The  further  stages  of  rough  entry  are  shown  in  Plate  6. 
Unfortunately  only  the  first  photograph  shows  the  sphere  in  the 
impact  phase  and  here  the  splash  cannot  be  seen  very  clearly. 
However,  the  sheath  can  be  seen  as  a  light  band  stirrounding  the 
sphere  above  the  surface. 

For  the  usual  projectile  entering  water  the  speed  does 
not  exceed  five  or  ten  percent  of  sound  speed  in  water.  It  is 
therefore  reasonable  to  neglect  compressibility  effects.  For 
the  sake  of  general  interest,  we  include  Plate  7,  which  shows 
a  sphere  entering  at  roughly  half  the  speed  of  sound  in  water. 
Here  the  effect  of  the  compressibility  is  of  course  inportant. 

Section  3.  Virtual  Mass 

The  impact  drag  coefficient  Cp  will  be  obtained  in  terms 
of  the  nvirtual  mass”  M  of  the  fluid.  Let  U  be  the  downward 


velocity  of  the  sphere  at  any  instant  t,  and  let  P  be  the 
resultant ' upward  impact  force  exerted  on  the  sphere.  Then 
the  sphere  exerts  a  downward  force  P  on  the  fluid  and  im¬ 
parts  momentum  to  the  fluid*  The  "virtual  mass"  M  of  the 
fluid  at  a  given  instant  t  is  defined  by 


(3*1) 


3£{MU) 


./pdt 


Let  M0  be  the  mass  of  the  incoming  projectile*  Neg¬ 
lecting  gravity  compared  to  the  impact  force  Pp  we  have 

(3.2)  -  P  =  ^(M0U)  . 


Adding  (3*1)  and  (3*2)  we  obtain  ^  [(M  +  MQ)  U]  =*  o  and  so 
integrating. 


(3.3)  (M  +  MjU  =  MU  or  U  =  - 2 - 

1  +  M/*0 


where  UQ  is  the  velocity  of  the  incoming  sphere  at  the  initial 
instant  of  contact  with  the  fluid* 

Let  B  =  B(t)  be  the  distance  at  time  t  of  the  lowest 
point  of  the  sphere  from  the  initial  surface  level;  the  velocity 
U  of  the  sphere  is  then  B  «  •  By  differentiating,  (3.3)  with 

*  Thus  MU  is  the  vertical  impulse  contributed  to  the  fluid  by 
the  moving  sphere.  It  is  also  possible  to  introduce  a  virtual 
mass  based  on  the  kinetic  energy  of  the  fluid,  but  for  the  case 
of  impact  it  differs  from  the  virtual  mass  defined  above.  This 
Is  to  be  expected  from  general  considerations  of  Inelastic  im¬ 
pact.  We  shall  use  throughout  this  report  the  definition 
(3.1)  of  the  virtual  mass,  based  on  momentum;  Justification 
of  this  choice  Is  given  in  the  Appendix. 


respect  to  t  and  using  (3*2),  it  follows  that 


( 3*4 ) 


P  =  jj2  = 

1  +  ® 

Mo 


dM/dB 

¥ 


(l  + 


in  \  i 

V 


Let  A  denote  the  radius  of  the  incoming  sphere,  and  let  p  be 
the  density  of  water.  Introducing  the  dimensionless  quantities 

(3*5)  m  =  -x  f  b  =  ~  , 

A 


we  may  write 

(3.6)  P  =  iyoU2  .  tr  A2  .  C 

^  O  Jr 


where 


For  incoming  projectiles  of  largo  mass  compared  to 
the  mass  of  the  water  displaced  by  the  hemispherical  nose, 
may  he  neglected  compared  to  1.  Then  the  inpact  drag 
c8ef ficient  Cp  is  merely -jg.  Also,  from  (3*5),  the  change 
in  velocity  U  in  the  impact  stage  is  small* 

However,  if  M0  is  not  large  compared  to  M,  the  denom¬ 
inator  in  (3*2)  must  he  retained*  This  is  the  case  in  the 
experiments  which  will  he  cited  in  Section  9,  and  it  is  con¬ 
venient  to  put  the  denominator  in  a  slightly  different  form. 
Let  <r  he  the  specific  gravity  of  the  incoming  projectile 
considered  as  though  the  whole  mass  were  concentrated  into  a 
sphere  (if  the  incoming  projectile  is  not  already  a  sphere), 
1  *  e* , 

(3.8)  <T  =  — Hh —  . 


Using  (3.5), 
(3.9) 


equation  (3.7)  may  he  written 

dm 

C  - - -  . 


The  determination  of  the  impact-drag  coefficient  Cp  requires 
merely  the  determination  of  the  dlmons ionless  virtual  mass  m* 

Section  4*  Hydrodynamic  Formulation  of  the  Problem. 

W©  shall  neglect  the  effects  of  compressibility  and  of 
the  viscosity  of  the  fluid.  Since  the  fluid  is  at  rest  before 
the  sphere  strikes  it,  the  flow  of  the  fluid  at  each  Instant 
of  time  t  is  irrotational  and  can  be  described  by  a  velocity 
potential  V  •  This  is  a  function  t/(x,  y,  z,  t)  defined  over 
the  region  of  space  occupied  by  the  fluid  such  that 

-  grad  xy 

is  the  velocity  of  the  fluid  at  the  instant  t  and  the  position 


x,  j,  z.  The  velocity  potential,  on  account  of  the  incompress 
ibillty  of  the  fluid,  satisfies  the  potential  equation 


^XX  +  Vyy  +  Vzz  =  0  , 

where  subscripts  mean  differentiation.  The  potential  is 
normalized  by  setting  Y  =  0  at  e»  ),  The  pressure  in  the  fluid 
is  connected  with  the  potential  function  by  the  Bernoulli 
equation 

(4*1)  £  .f  i  (gradl J/)^  -  =  0 

where  p  is  the  pressure  in  excess  of  hydrostatic  pressure* 


i 

|  z  -ax  /s 
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To  obtain  the  Impact  force  exerted  on  the  sphere, 
we  shall  use  the  principle  of  conservation  of  momentum. 
The  total  momentum  of  the  fluid  is  directed  vertically 
downward  because  of  the  rotational  symmetry  and  is  equal 
to 


(4*2)  <3bcdy&z  as  fjjv dxdy  -pjjip'dx&y  • 

w  *  s  //«, 


Here  W  is  the  wetted  portion  of  the  sphere,  S  is  the  free 
surface,  and  the  last  term  in  (4# 2)  means  the  limit  of  the 
integral  JJy dxdy  extended  over  the  submerged  portion  H  of  a 
large  spherical  surface  the  radius  of  which  tends  to  infinity. 

The  external  forces  acting  on  the  fluid  contained  in 
the  membrane  H  are:  (i)  the  downward  force  of  magnitude  P 
exerted  by  the  sphere  and  (ii)  the  upward  force  exerted  on  H 
by  the  pressure  of  the  fluid  external  to  it.  (The  resultant 
force  due  to  the  constant  atmospheric,  pressure  is  zero.  The 
only  effect  of  gravity  is  the  buoyancy  which  is  negligible 
compared  to  the  intact  2*orc:e  P.  )  Tha  force  (li)  is  equal  to 
th®  time  rate  of  change  of  the  last  term  in  (4.2).  For,  by 
Bernoulli's  law  (4.1),  the  force  (ii)  differs  from  the  time 
rate  of  change  of  the  last  term  by  a  term  of  the  form 


//(grad  >50 2 

H 


dxdy 


which  approaches  zero  as  the  radius  of  H  becomes  infinite. 
(For,  at  co  ,  y  behaves  at  worst  like  -i  , therefore  (grady)2 
like  — ^  and  the  whole  integral  like—g)  •  Likewise  the 
momentum  flux  through  the  surface  H  consists  of  integrals  of 
squares  of  velocities,  and  these  approach  zero  as  the  radius 
of  H  becomes  infinite.  Since  the  resultant  of  all  the  ex¬ 
ternal  forces  is  equal  to  the  time  rate  of  change  of  the 
momentum  (4.2),  wo  therefore  have 


* 


(4.3) 


p  =  (f/fvteccy) 
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The  formula  (4*3)  can  also  he  Inferred  directly  from 
the  known  physical  interpretation  of  the  potential  function 
y  ,  namely  that  py  is  the  impulsive  pressure  required  to 
produce  the  flow* 

It  is  convenient  to  exhibit  separately  the  effect  of 
the  velocity  B  of  the  sphere  by  setting 

(4.4)  7p  z»  t)  =  B<|)(x,  y,  z,  t)  , 

so  that  is  the  instantaneous  potential  function  correspond* 
ing  to  the  parameter  B  and  to  unity  as  velocity  of  the  sphere* 
Then  (4.3)  can  be  written 

p  -&  <®> 

where 

(4.5)  U  -  • 

W+S 

This  formula  expresses  the  virtual  mass  of  the  fluid  in  terms 
of  the  potential  function  $  • 

Section  5*  The  First  Approximation  for  a  Sphere 

Consider  the  flow  produced  in  the  fluid  as  a  result 
of  the  entry  of  the  sphere*  As  the  sphere  penetrates  the 
surface,  the  free  surface  no  longer  remains  plane  but  rises 
slightly  and  forms  spurts  running  up  the  sides  of  the  sphere. 
Although  the  exact  shape  of  the  surface  is  difficult  to  deter¬ 
mine  theoretically,  some  estimates  will  be  presented  in  Sec¬ 
tion  6*  Por  the  purpose  of  computing  the  virtual  mass  and  the 
impact  force,  we  shall  make  a  plausible  first  approximation 
in  the  present  section  by  disregarding  the  rice  of  the  surface* 
More  exactly,  we  shall  neglect  the  rise  of  the  surface  and 
also  the  squares  of  velocities  on  the  surface* 


To  determine  the  potential  function  y  (or  <f>)t  it 
is  necessary  to  know  boundary  condition  for  y  •  The  exact 
boundary  condition  on  the  free  surface  S  of  the  fluid  is 

P  =  0 

(pressure  is  atmospheric),  which  by  Bernoulli fs  equation 
(4*1)  can  be  expressed  in  terms  of  the  potential  function 

Y 

(5.1)  |T  s  1  (grady)2  . 

Since  this  exact  boundary  condition  is  non-linear  and  diffi¬ 
cult  to  apply,  we  neglect  the  rise  of  the  surface  S  and 
squares  of  velocities  on  S  as  a  first  approximation.  Thus, 

O 

neglecting  (gradV')  in  (5.1)  we  have 

§£=  0 

or,  since  the  free  surface  S  is  assumed  to  remain  at  its 
original  level, 

( 5® 2 )  y  =  0  on  the  initial  plane  z  =  0  (since  y^  0  at  »  )• 

The  boundary  condition  (5*2)  states  that  the  impulsive 
pressure  f>y  is  zero  on  the  surface.  This  condition  would  be 
correct  if  the  impact  took  place  completely  impulsively i 
Actually  the  time  interval  is  very  short,  and  so  (5.2)  is  a 
good  first  approximation. 

It  is  more  convenient  to  use  the  potential  function 
$  for  unit  velocity,  defined  by  (4«4).  The  potential 
function  $  (depending  onB)  satisfies  the  following  conditions 


(a)  On  the  sphere,  the  negative  of  the  normal  derivative 
of  y  ss  B(f)  is  equal  to  the  component  of  the  velocity  of 
the  sphere  along  the  normal,  or  ' 

-  =  cos  a  • 

(b)  On  the  initial  water  level, 

(f)  =  0  » 

There  Ls  a  unique  potential  function  satisfying  these  con¬ 
ditions. 

The  flow  arising  from  the  above  boundary  conditions 
has  been  constructed  explicitly  In  the  preceding  report  [15]. 
Because  of  condition  (b),  the  potential  flow  can  be  extended 
to  the  infinite  space  by  reflection  on  tho  surface  and  the 
flow  at  each  instant  corresponds  to  the  steady  flow  about  a 
symmetrical  lens  formed  by  the  Intersection  of  two  equal 
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spheres.  A  photograph  from  Worthington  [19]  is  shown  "below. 
This  photograph  is  a  view  of  the  entering  sphere  from  slightly 


below  the  surface.  The  surface  reflection  creates  the  upper 
side  of  the  lens.  (Of  course,  the  size  and  shape  of  the  lens 
varies  with  time.) 

The  preceding  report  Ll5]  supplies  not  only  the  potential 
function  <t>  but  also  the  corresponding  virtual  mass,  which  we 
shall  henceforth  denote  by  the  subscript  L  referr lag  to  a 
lens.  The  results  of  [15]  which  we  use  here  may  be  r  uumarized 
as  follows: 

(i)  The  dimensionless  virtual  ma&o  *=  — -L —  is  a  function 


1-..3 
t?*  p  A 


of  b  alone  (where  b  =  B/A),  and  will  be  written  as  M^(b).  Its 

The  impact-drag  coefficient 

The 


dependence  on  b  is  drawn  in  Graph  3. 
dror(b) 

— gr —  depends  on  b  in  the  manner  indicated  in  Graph  2. 


maximum  of  is  .95  and  It  occurs  for  b  =  .24 • 

When  b  is  small,  the  spherical  segment  beneath  the  initial 

level  of  the  fluid  is  nearly  a  disc*  By  a  known  formula  for  the 

virtual  mass  M  of  a  disc  (see  [8],  p.  130)  we  therefore  have 

,,4  3 

M  =  7/oc 
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for  B  small*  where  c  Is  the  radius  of  the  circle  of  Inter¬ 
section  of  the  sphere  with  the  original  level. 


This  gives 


»L<b) 


M 


5V1' 


=  ^  (2b-b2)3/2 


or 


(5.3) 


m^b)  *  23/2  b3/2  =  2.40  b3/2 


for  small  b.  The  more  exact  lens  flow  gives 

m^b)  =  2.40  b3/2  -  1.15  b2  +  ... 

=  3.60  b1/2  -  2.30  b  +  ... 


(5.4) 


db 


(11)  The  variation  of  the  virtual  mass  with  the  depth  of 
penetration  can  be  better  visualized  by  setting 


(5.5) 


Ml  -  kjj{°  ©  t 


where  the  dimensionless  quantity  depends  on  b.  The  coefficient 
^L  *3  drawn  in  Graph  5.  When  b  Is  small  the  spherical  segment 
is  nearly  a  disc  so  that 


3  aa  b 


0  . 


When  b  =  1,  the  spherical  segment  is  a  hemisphere  and  It  Is  then 
well  known  that  |  ( see  [ 8] ,  p.  116).  The  quantity  kL  de¬ 

creases  rapidly  from  the  value  4  at  b  =  0  and  then  levels  off. 

oven  swinging  upward  slightly  near  b  =  1.  This  behavior  of  kT 

jl» 

is  plausible  on  Intuitive  grounds.  (The  final  coefficient  k, 
which  includes  the  various  corrections  made  in  this  report,  is 
also  drawn  in  Graph  5. ) 


It  is  interesting  to  note  that  for  the  case  of  a  cone 
Instead  of  a  sphere  the  quantity  k  is  a  constant  independent 
of  the  depth  of  penetration  of  the  cone.  This  results  from 
the  fact  that  for  a  cone  the  flow  patterns  at  different  times 
are  geometrically  similar. 

To  obtain  a  second  approximation  to  the  virtual  mass  and 
to  the  impact  force,  it  is  necessary  to  estimate  the  effects  of 
the  motion  of  the  surface.  On  intuitive  grounds,  corrections 
due  to  the  motion  of  the  surface  can  be  divided  into  two  parts. 

On  the  one  hand  the  water  rises  and  wets  a  larger  portion  of 
the  sphere,  the  sphere  now  exerting  pressure  on  the  fluid  over 
a  larger  area  and  causing  an  Increase  in  the  total  Impact  force. 
This  will  be  called  the  "wetting"  correction.  On  the  other  hand, 
the  free  surface  of  the  fluid  rises  and  relaxes  somewhat  the 
restraint  imposed  on  the  sphere.  This  will  be  called  the  "stir- 
face"  correction. 

These  two  corrections  appear  in  formula  (4.5)  when  written 


as 

(5.6) 


M 


o JJ  $)  dxdy  +  ^ j'J  $  dxdy 


W  S 


The  second  integral  is  not  zero  since  the  condition  (5.1)  re¬ 
places  $  =  0.  (In  fact,  $  will  turn  out  to  be  negative  on  S. ) 


Section  6.  The  Wetting  Correction 

The  greatest  disturbance  of  the  surface  of  the  fluid  occurs 
alongside  the  sphere.  Here  a  thin  sheath  of  the  fluid  rises  more 
rapidly  than  the  rest,  and  may  detach  from  the  sphere  (see  Plate  4). 
This  sheath  because  of  Its  thinness  contributes  very  little  to 
the  impact  force  (the  pressure  of  the  fluid  in  the  sheath  is  practi¬ 
cally  atmospheric).  The  sphere  exerts  pressure  on  the  fluid  in 
the  main  over  the  portion  of  the  sphere  obtained  by  cutting  off 
the  thin  sheath  from  the  thick  "base"  of  the  risen  fluid.  This 
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portion  of  the  sphere  is  indicated  by  the  points  QQ  in  the 
diagram,  and  the  separation  of  thin  sheath  from  the  thick 
base  is  indicated  by  the  dotted  lines  QQ*.  The  flow  of  the 
fluid  is  due  to  the  portion  QQ  of  the  sphere,  and  is  approxi¬ 
mately  the  flow  about  a  lens  consisting  of  the  spherical  seg¬ 
ment  QQ  and  a  horizontal  plane  free  surface  QQ’oo  •  (This 
may  in  fact  be  taken  as  the  matnematical  definition  of  the 

•g. 

point  of  demarcation  Q  between  sheath  and  base.)  In  other 
words,  the  "base0  below  QQ’  is  so  thick  that  its  influence  on 
the  sphere  is  practically  the  3&me  as  if  it  extended  to  oo  , 

It  will  be  shown  later  in  this  section  that  the  above  discussion 
is  exactly  verified,  and  the  definition  of  Q  justified,  in  the 
asymptotic  case  as  B  — >  0  • 


*  For  the  method  used  In  the  practical  determination  of  this 
point  of  demarcation  Q  from  photographic  data,  see  Plate  4 
and  the  accompanying  discussion  on  page  2-3. 
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Henceforth  it  will  be  the  spherical  segment  QQ  which 
will  be  called  the  "wetted"  portion  W,  or  the  "effective" 
portion,  of  the  sphere.  Let  be  the  depth  of  this  effective 
portion,  and  the  radius  of  its  circular  base  (see  diagram). 
The  quantity  will  be  called  the  "effective  depth"  of  penetra¬ 
tion  of  the  sphere*  We  shall  also  set 


b  _  B1 
bi  ‘X  • 


The  flow  of  the  fluid  is  due  to  the  effective  portion  of  the 
sphere  and  is  approximately  due  to  a  lens  with  dimensions  B^, 
in  place  of  B,o. 

The  dimensionless  virtual  mass  m  =  — & —  is  given  by 


2  P  A 

a^(b^)  of  the  lens  flow*  Thus  the  dimensionless  virtual  mass 

and  the  drag  coefficient  ~  are 

db 


(b.l) 


m 


=  “L<bl) 


dm_dtaL(b1) 
db  db-, 


dm-Cb,) 

where  m^(b^)  and  ■  — -  —  tire  determined  from  the  lens  flow 
with  the  argument  b,. 

It  remains  then  to  determine  b^  as  a  function  of  b.  Set 
(6.2)  b1  =  b  .  w(b ) 


(or  equivalently,  B^  =  B*w(b)  )  where  w(b)  is  a  function  of 
b  =  B/A.  The  quantity  w(b)  will  be  called  the  "wetting  factor." 
In  this  section  we  shall  show  theoretically  that  initially, 
for  b  =  0,  the  wetting  factor  w(0)  has  the  value  3/2.  The  exact 
dependence  of  w(b)  on  b  la  extremely  difficult  to  determine 
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theoretically,  hut  for  small  values  of  b  the  theory  shows 

*T  . 

that  w(b)  behaves  like  ^  -  constant  ./b^  Consequently, 

the  variation  of  w(b)  with  b  was  determined  from  experimental 

measurements* 

A  blackened  steel  sphere  of  diameter  2*86  cm  was  dropped 
into  a  tank  of  water  from  a  height  of  16*2  cm  and  flash  photo- 

■jj- 

graphs  were  taken  with  the  sphere  at  various  stages  of  immersion. 
The  photographs  in  Plates  2  and  3  are  a  selection  showing  the 
sequences  of  events.  Plate  4  shows  an  enlarged  view  of  the 
'  splash  and  illustrates  how  the  measurements  were  made*  The 
straight  line  represents  the  Initial  surface  level,  the  posi¬ 
tion  of  which  was  determined  by  a  pair  of  inclined  needles  and 
their  reflections,  as  indicated  in  the  photograph.  The  boundary 
of  the  sphere  has  been  completed  by  a  white  line.  The  exact 
line  of  demarcation  between  sheath  and  base  is  somewhat  arbitra¬ 
ry  but  estimates  can  be  made  with  a  fair  degree  of  accuracy  and 
consistency*  The  upper  pair  of  short  horizontal  lines  shows  our 
choice  of  demarcation  between  sheath  and  base.  Below  this  line 
the  risen  water  is  thick  while  above  it  the  water  thins  out  very 
rapidly.  For  this  particular  photograph,  the  values  of  b  and 
w(b)  are: 

b  =  .27 

w(b )  =  1.30 

Similar  measurements  were  made  on  a  number  of  enlarged 
photographs  showing  various  depths  of  penetration.  The  quantity 
w(b)  was  then  plotted  against  /b  with  the  results  indicated  in 
Graph  8(a).  Also  clotted  in  Graph  8(a)  is  an  experimental  point 
obtained  by  Watanabe  while  measuring  impact  forces.  This  last 
point  is  discussed  in  Section  9. 


#  The  experiments  were  performed  at  the  AMG-NYU  Laboratory  by 
M«  Shamos. 


The  experimental  points  shown  in  Graph  8{a)  fall  on 
a  reasonably  well-defined  straight  line  passing  through  the 
theoretical  value  w(0)  =  3/2.  The  equation  of  tills  line  is 

(6» 3)  w(b)  =  1*5  —  *4  /b  . 

In  Graph  8(b),  w(b)  is  drawn  as  a  function  of  b.  Actually, 
in  Graph  8(a),  the  experimental  points  are  all  slightly  above 
the  straight  line  (6*3)  near  /IT  =1*  A  more  accurate  curve 
fitting  the  data  is  the  parabola 

w(b)  =  1.5  -  .4/b  +  .065  b  . 

But  this  will  make  a  noticeable  difference  only  near  b  =  1,  and 
we  shall  rather  use  the  straight  line  (6.3). 

It  remains  to  derive  the  theoretical  value  w(0)  =  3/2  . 
We  shall  use  a  method  similar  to  one  applied  by  Wagner  [16] 
for  the  case  where  the  entering  body  is  a  nearly  flat  wedge 
(two  dimensional) *  In  this  case,  the  wetting  factor  turns  out 
to  be  £  =  1.57.  A  corresponding  calculation  for  a  very  flat 
cone  gives  a  wetting  factor  —  =  1.27. 

Tf 

Consider  the  wetting  correction  for  small  b,  when  the 
sphere  has  entered  only  a  small  distance  Into  the  fluid.  The 
flow  of  the  fluid  is  approximately  due  to  the  spherical  segment 
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QQ.  Since  B  and  sire  supposed  small  the  flow  due  to  this 
spherical  segment  can  be  approximated  by  that  due  to  a  cir¬ 
cular  disc  of  radius  c^.  These  approximations  are  accurate 
in  the  limiting  case  as  B  — >  0,  as  can  be  seen  by  -the  follow¬ 
ing  argument*  Perform  a  similarity  transformation  with  magni¬ 
fication  factor  -i  .  The  limiting  situation  as  B  -►  0  is  indi  - 

c 

cated  by  the  diagram  below.  The  effective  portion  of  the  sphere 


in  contact  with  the  water  becomes  the  circular  disc  QQ,  while 
the  free  surface  becomes  the  plane  Q  oo  . 

The  theory  of  the  flow  due  to  a  circular  disc  is  classi¬ 
cal  (  [  8]  *  pp.  130,  135,  or  [10],  pp.  456,  457).  The  potential 
function  yf  describing  the  flow  is  given  in  [8],  p.  135.  In 
particular  the  value  of  y  on  the  lower  surface  of  the  disc  is 


(6.4) 


t 


where  r  is  the  distance  from  the  center  and  c^ 
of  the  disc.  The  value  of  the  upward  velocity 
free  surface  of  the  fluid,  level  with  the  disc. 


is  the  radius 

-  on  the 

dz 

is  (  [ 8] ,p.  130): 


-  sr 


3z 


(6.5) 
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The  motion  of  the  surface  of  the  fluid  can  now  be 
followed.  At  any  point  of  the  free  surface  at  distance 
r  from  the  axis  of  symmetry,  let  7\  (r)  be  the  height  of  the 
free  surface  at  the  time  t  above  its  original  level  (the 
dependence  of  *j(r)  on  the  time  t  is  understood*  Assume  that 
the  upward  velocity  of  the  surface  at  the  instant  t  is  that 
Induced  by  the  disc  of  radius  c^*  Prom  (6*5),  we  have 


The  quantity  c^  is  a  function  of  B^,  the  relation  between 
them  depending  on  the  shape  of  the  incoming  object.  For  a 
sphere  of  radius  A  we  have 

(6.8)  of  »  2ABX  -  B^  . 


o 


(The  upper  limit  c^  of  Integration  should  be  smaller  than 
r.  The  fluid, particle  initially  at  a  distance  r  rises  until 
it  hits  the  sphere,  after  which  it  enters  into  the  thin  sheath)* 
In  this  integration  r  is  kept  fixed*  Set 

(6.11)  -1  »  t 

V 


where  £  -  1.  Then  (6*10)  becomes 

.2 


*?(r)  - 


_  r 


Aw(0) 


F(S)4d£ 


or,  using  (6.11)  and  (6*9), 
(6.12)  -n(r)  =  B  . 


I  d|. 


Equation  (6*12)  represents  the  equation  of  the  surface 
as  long  as  r  -  c-^.  In  particular,  at  the  position  r  =  c^,  or 
£  *=  1,  the  height  of  the  surface  is 


B1  -  B  =  B  .  Cw(0)  -  1]  . 

Equation  (6.12)  yields 

(6.13)  w(0)  -  1  =  2f  nl  1  |  d  \  . 

Substitution  of  (6*7)  into  (6*13)  gives  the  value  of 
the  wetting  factor  w(0)  at  b  =  0: 

(6.14)  w(0)  =  |  . 


Evaluation  of  (6.12)  gives  the  equation  of  the  surface  ©£~the- 
water 


V  (r)  =  B  •  i  [  {—■ 2  -  2)  arc  sin  ^  -  -  y4  -  £2]  . 


(6*15) 
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Equation  (5*14)  is  plotted  below  for  the  case  when  b  -  «01. 
(Of  course,  this  method  cannot  be  expected  to  produce  the 
thin  sheath* ) 


Section  7.  The  Free  Surface  Correction 

The  free  surface  correction  arises  from  the  tena 

(7.i)  fjfyoxitj 

s 

occurring  in  (4*3)  and  the  fact  that  w©  shall  no  longer  neglect 
squares  of  velocities  on  3*  The  exact  boundary  condition  is 
then 


(7.2) 


ft  =  2  (?raa''*'>2 


in  place  of  y  -  0  (see  (5*1)  and  the  accompanying  discussion). 
Suppose  we  follow  a  particle  of  fluid  on  the  free  surface 
and  determine  the  rate  at  which  the  value  of  V  changes  for 
that  particle.  Thl °  requires  use  of  the  "particle  derivative" 

where 

Dt 

=  ■!?  “  grady  .  grad  y  . 
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In  view  of  (7.2),  we  have 

(7.3)  jjjf  =-  |  (grad^)2 

on  the  free  surface  S.  The  value  of  Y  on  S  thus  constantly 
decreases  and  is  negative  at  any  instant  after  the  initial 
contact  of  the  sphere  with  the  fluid.  This  shows  that  the 
surface  correction  (7.1)  is  negative. 

The  correction  to  the  impact  force  P  arises  from  the 
time  derivative  of  (7.1).  Let  the  height  of  the  surface 
above  the  initial  level  be 

V  =  V  (r,t) . 

Denote  the  surface  correction  to  the  Impact  force  by  P5  ,  so 
that  from  (4.3), 

(7.4)  Ps  =  tody  =/“//{§£  +  ■n  If)  tody  . 

5  S 

(There  is  no  contribution  from  the  changing  domain  of  Integra¬ 
tion  since  this  is  cancelled  by  the  corresponding  term  from 

( /°//l^dxdy)  •)  Substitution  of  (7.2)  gives 
at:  w 

(7.5)  Ps  */5jJ'||(grad'^)2  +  V  |^|  dxdy  • 

s 

Now  assume  as  a  reasonable  approximation  that  the 

velocity  on  the  surface  S  is  that  induced  by  the  lens  flow, 

used  as  a  first  approximation  in  Section  5.  The  velocity  of 

the  surface  is  then  completely  vertical,  so  that  nj  =  -  ^ 

o  z 

and  (7.5)  becomes 

(7 ..6)  Ps  =  -  |/»//’(||)2  tody  =  -  |  B^ff)2  tody  . 
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The  velocity  -  is  infinite  at  the  corner  of  the  lens, 

dz, 

hut  the  integral  (7*6)  nevertheless  converges*  However, 
the  integral  (7*6)  is  a  bit  cumbersome  to  evaluate  for  a 
lens  flow;  we  shall  instead  approximate  the  flow  at  each 
instant  by  that  induced  by  an  oblate  spheroid  circumscrib¬ 
ing  the  lens*  (In  the  previous  report  [5],  this  oblate 
spheroid  was  shown  to  be  a  reasonable  approximation  to  the 
lens* ) 

The  flow  about  an  ellopsoid  is  well  known  (see  [10], 
pp*  456-7)  and  [8],  p*  132  ff*  )• 
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7* 


where 


oos  'i 


Kg  1b  drawn  as  a  function  of  b  =  B/A  In  Graph  7.  This  correc¬ 
tion  Is  not  to  be  considered  as  accurate  for  very  small  b,  but 
should  be  a  reasonable  estimate  over  mos^  of  the  range  In  which 


we  are  Interested* 


Section  8*  The  Second  Approximation  1 

The  results  of  Sections  5,6,7  will  now  be  collected  to¬ 
gether*  Set 

\  “/’//Wxd i  > 


so  that  is  the  contribution  to  the  virtual  mass  (4*5)  due 
to  the  wetted  region  W  alone.  Then  equation  (4*3)  can  be 


written  in  the  form 


P  =  ~ 


<V>  +  p, 


where  P  is  given  by  (7*4),  or 


»  diL  .  2 
?  "  V  +  -3?  B  + 


Using  MqB  =  -  P  and  (7*7),  we  have 


l.iZ  _  .zf*3. 


m  +  3=  )  =  i.B  .  tr  A*  -  K 

Un  2  '  dh  J 


=  ipU2  .  tr  A2  . 


P  =  ±  pTT 
2  r  o 


»  - 


(i  +  ?S)(1  +  M  f 


(8.1) 
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M  \3 


In  this  formula,  the  denominator  can  be  approximated  by  (1+  ~  )  * 

Mo 

Comparing  fonmila  (8*1),  the  denominator  being  replaced 
by  (1  with  formulas  (3*6),  (3*7),  *9  see  that 


(8.2) 


^  a  -J 

db  db 


where;  m  is  the  final  estimate  of  the  virtual  mass,  and 


=  W  • 


Thna  m  may  be  obtained  by  integration,  with  the  result 


(8.8) 


“L(bl> 


K  db 

s 


This  is  drawn  in  Graph  3  and  compared  with  the  first  approxima¬ 
tion  Bj^b)  arising  from  the  lens  flow. 

The  Impact  drag  coefficient  is  calculated  from  formula 

(3.7)  using  the  final  estimate  of  su  The  numerator  ^  can 

db 

be  written,  by  (8.2),  in  the  form 

(8.4)  a  -  .  Si  -  K  . 


dm- (b. ) 

The  quantity -  may  be  read  from  Graph  2,  which  is  a  repeti- 

db. 

1  db. 

of  the  result  obtained  in  [15]j  — is  obtained  from  formulas 

(jb  -a 

(6*2),  (6*3):  and  K.  is  read  from  Graph  7.  The  result  for  ~ 

8  db 

is  drawn  In  Graph  1  (the  case  <r  =  00  ).  This  is  exactly  the 

impact  drag  coefficient  (V  in  case  the  mas3  Mq  of  the  Incoming 

projectile  is  large  compared  to  the  mass  of  the  water  displaced 

by  the  hemispherical  nose. 
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The  impact  drag  coefficient  Cp  in  the  general  case 

of  finite  s'  may  be  calculated  from  formula  (3*9)  and  tbe 

graphs  of  m  and  •==•  The  results  corresponding  to  c  =  °°  , 

db 

1,  #5,  *25  are  shown  in  Graph  1* 

We  observe  that  the  maximum  value  of  Cp  is  about  1 
(unless  <r  is  wary  small),  and  that  this  maximum  occurs  for 
b  =  *15  approximately* 

In  the  neighborhood  of  b  =  0  the  expansion  of  Cp  can 
be  obtained  from  (8*4),  using  (5*3)  or  (5*4).  The  result  is 

(8.5)  Cp  *  6.62  b1/2  -  a2b  +  ... 

where  a2  is  a  coefficient  which  is  not  accurately  given  by 
(8.4)  and  its  determination  is  difficult.  The  coefficient 
6.62  Is  i  .  33/2. 

tf 

In  Graph  2,  the  final  Cp  is  compared  with  the  Cp  ob¬ 
tained  from  the  first  approximation  of  the  preceding  report [15]. 


Section  9.  Comparison  with  Experiment 


Careful  experiments  on  the  impact  of  cones  and  spheres 
have  been  conducted  by  S.  Watanabe  (see  [17).  He  used  a  thin 
spherical  segment,  with  dimensions  as  indicated  in  the  diagram, 
on  top  of  which  were  placed  a  weight  and  piezoelectric  gauges 
connected  to  an  oscillograph.  The  system  was  dropped  from 
heights  of  20  cm  to  50  cm  from  the  water  surface,  and  the 
oscillogram  record  gave  directly  the  impact  force  as  a  function 
of  the  time.  The  Watanabe  results  arc  exhibited  In  Graph  8(a), 

\ 

\ 


<r 


h 

Maas  Mq 


=  ..84  cm.  \ 

=1.645  kgm.  ' 

^  /S'  c  rn . 

=  .116 


T  h  =  cm  . 
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taken  from  his  paper*  the  vertloal  axis  R  "being  the  Impact 
foroe  measured  In  kilograms  of  force  and.  H  being  the  height 
from  vhloh  the  spherical  segment  was  dropped* 

These  graphs  are  converted  Into  graphs  of  Cp  as  follows: 


jr/^o  *1,a2  §/*•' 


WyoA2H 


b  =  depth  of  Immersion 

A  ' 

where  the  radius  A  of  the  sphere  Is  15  cm*  Tho  results  are 
exhibited  In  Graph  8(b)* 

Theoretically  one  should  expect  the  formula  (3*9)  for 
Cp  to  be  valid*  Expanding  m  and  <-|S  In  powers  of  b* 


(0.1) 


»  *  |  ,lbV8  -  4?  b2  +  ... 


tin 

•IM  SS  fl 

db 


.  a  b  + 


where  a1#  ag  are  constants  (from  (8.3)  the  theoretical  estimate 
of  a^^  is  6*62.  The  denominator  (1  +  ~)3  «  l  +  2®  near  ^  =  o 
so  that  ec  8<r 


(9.2) 


a-,bV2  .  a  b 

C  «,  - - - - 

1  +  3  aib3//2  _  _9_  a2b< 
4  <r  16  <r 


1  =  ^/8  [i-fk] 


°p  16 


3  • 
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Sotting 


(9.3) 


wo  have 

(9.4)  V  *  a2*  “  1  - 

The  quantities  y  and  x  are  plotted  in  Graph  9(a) 
using  Watanaoe^  experimental  points.  They  show  an  excellent 
agreement  with  the  linear  relation  (9.4),  and  the  best  straight 
lines  are  likewise  drawn  in  Graph  9(a),  the  resulting  values 
of  a^  and  ag  being: 


H(  cm) 

20 

30 

40 

50 

al 

6.96 

6.92 

7.17 

7.13 

a2 

23.5 

25.9 

26.6 

27.3 

The  values  of  a^  show  good  agreement  with  the  expected  theoreti¬ 
cal  value  6*62. 

Selecting  in  equation  (9.2) 

(9.5)  a1  «  6.62  , 

and  dropping  the  denominator  completely,  we  have 

6.62  b 1^2  -  e 
- R 


(9.6) 


b 
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6.62  b1/2  -  c 

The  quantity  - - -P  is  plotted  against  b  In  Graph  9(b), 

the  values  of  b  less  than  .01  being  excluded  since  any  experi¬ 
mental  error  in  c^  will  then  produce  a  large  error  in  ag. 

The  result  shows  that  the  quantity  (9.6)  is  accurately 
constant  and  its  best  value  from  the  experimental  data  is  23.5. 
If  the  denominator  in  (9.2)  is  included,  a  correction  to  the 
value  of  &2  is  obtained  which  amounts  to  -  1.5.  Thus,  the  ex¬ 
perimental  value  of  &2  is 

t 

(9.7]  a2  =  22  . 

Thus  the  theoretical  result  (9.5)  Is  verified  completely 

3 

by  Watnnabe*s  experiment.  In  particular,  the  value  —  for  the 
wetting  correction  is  verified. 

Another  more  direct  check  of  the  wetting  correction  of 

X 

g  is  provided  by  Watanabe  In  his  observation  of  the  depth  of 
immersion  at  which  the  resistance  curve  dropped  sharply.  The 
depth,  in  Table  IV  of  Watanabe‘s  paper  [17],  is  given  as  .585cm 
on  the  average.  This  is  clearly  due  to  the  water  overflowing 
the  spherical  segment.  The  wetting  factor  is  therefore  approxi¬ 
mately 

w(b )  =  -A§±  a  1.44 

.585 


which  is  slightly  less  than  3/2.  The  corresponding  value  of 
b  is 


b  ts 


.585 

15 


.039  . 


This  experimental  point  is  plotted  in  Graph  6(a). 

In  the  experiments  by  Blundell  [3],  accelerometers  with 
low  natural  frequencies  of  100-400  cycles/sec  and  low  damping 
were  used.  The  responses  of  the  instruments  were  oscillatory. 
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»o  that  only  the  general  shape  of  the  Impact  curve  could 
he  established.  The  dimensions  of  the  dropped  sphere  sire 
as  Indicated* 


(A)  ii0  * 

(B)  M0  = 

(C)  Mq  = 

(D)  M0  = 


351b,  UG  =  11.4  ft/sec,  <T  =  .23 

321b,  U0  =  11.4  ft/sec,  <T  =  .21 

321b,  UQ  =  8  ft/sec  ,  <T  s  .21 

501b,  UQ  =  8  ft /sec  ,  <T  =  .33 


The  experimental  curves,  obtained  from  Figures  1,  2  of  [3]are 
converted  into  C  curves  in  Graph  11(a),  and  the  corresponding 

Jr 

theoretical  curves  are  drawn  in  Graph  11(b).  An  idea  of  the 
experimental  error  can  be  obtained  by  comparing  (B)  and  (C), 
which  should  be  the  same,  and  by  noting  that  (D)  should  give 
the  highest  Cp  since  <r  is  the  largest.  The  agreement  with 
theory  is  satisfactory,  and  curve  (D)  even  agrees  numerically 
with  theory. 

Richardson  [13]  used  a  hemispherical  shell  with  the 
same  dimensions  as  Blundell,  but  with  a  mass  of  24  lb  (<r=  .159) 
dropped  from  a  height  of  2  ft.  He  used  an  instrument  which 
worked  on  the  changes  in  capacity  of  a  "breathing*  condenser. 

His  experimental  graph,  Figure  2b  of  [13],  is  converted  into 
the  dimensionless  form  of  a  Cp  curve  and  compared  with  the 
theoretical  curve  in  Graph  12.  The  agreement  is  satisfaorory, 
especially  on  the  position  where  the  maximum  Cp  is  attained. 

In  all  these  experiments,  the  value  of  <r  is  rather  small. 
It  is  desirable  to  perform  further  experiments  with  a  larger  <r  , 
uaing  accurate  instruments  (such  as  piezoelectric  gauges). 
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Appendix 

Section  10*  Comparison  of  Energy  and  Momentum 

In  Sections  3  and  4  the  force  of  intact  was  derived 
freon  momentum  considerations,  and  we  found  that  the  force 
on  the  entering  sphere  is  equal  to  the  time  rate  of  change 
of  the  resultant  impulsive  pressure  force  acting  od  the 
boundary  of  the  fluid  (formula  (4.3)  ).  We  may,  however, 
derive  the  force  of  ispaot  from  eiiergy  considerations. 

Let  £  M.TJ2  be  the  kinetic  energy  of  the  fluid  at  time 
t,  where  is  the  "energy*  virtual  mass.  By  the  principle 
of  conservation  of  energy,  we  have 


Differentiating  with  respect  to  t  and  using  (3.2),  we  obtain 


(10*2) 


P  « 


dM 

1  JE# 


!♦* 


1  £ 
2 - 


d  +  ^i)2. 

«o 


by  (10.1). 

Formulas  (10.1),  (10.2)  should  he  compared  with  the 
corresponding  formulas  (3*3),  (3.4)  based  on  momentum  con¬ 
siderations.  An  imnediate  comparison  of  (10.1)  and  (3.3) 
shows  that 

(10.3)  l+-~=(l+J)2or  Me=2M(l  + 

Ho  Ho  ° 
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where  M  Is  the  Momentum11  virtual  mass  which  we  have  been 
using  throughout*  The  some  relation  is  obtained  by  compar¬ 
ing  (10*2)  and  (3.4).  For  incoming  projectiles  of  large 
mass  compared  to  the  mass  of  the  water  displaced  by  the 
hemispherical  nose,  ^  J^/M0  is  small  compared  to  1,  and  we 
have  approximately 

(10.4)  II  =  2M  • 

y 

Thus,  for  impact,  MQ  Is  different  from  M. 

In-order  to  obtain  an  expression  for  Me,  we  observe 
that  the  total  kinetic  energy  of  the  fluid  is 

|  a.®2  =  taas-dz  -  -  IffjvU  33  -  63 

by  Green* s  theorem,  where  the  differentiation  is  with  respect 
to  the  normal  pointing  into  the  fluid.*  Thus  the  kinetic 
energy  Is  equal  to  the  work  done  by  the  Impulsive  forces  acting 
on  the  boundary  of  the  fluid*  Since  -  ^  =  TJ  coscx.  on  W  (see 
diagram  on  page  17),  we  see  that  -  dS  =  U  dxdy.  Substitu¬ 
ting  from  (4.4),  we  therefore  obtain  for  H  the  formula 

6 

(10.5)  *,=/>//♦  dxdy  +  />//*  i-  i  )  dS  . 

w  S 

Comparing  the  right  hand  sides  of  (10.5  and  (4.6), we 
see  that  they  differ  only  In  their  second  terms,  the  free  surface 
terms*  If  the  same  first  approximation  is  made  as  previously, 
that  is  $  =  0  on  the  surface,  we  would  obtain 

Me  =  fff*  t*6*  =  M 

- - - & _ _ 

*  It  is  easily  seen  (from  the  considerations  of  Section  4)  that 
the  integration  extended  over  a  large  spherical  membrane  tends 
to  zero  as  the  membrane  tends  to  infinity. 
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whereas  we  know  that  M.  =  2M  approximately*  The  comple te 
negleot  of  the  surface  terms  may  thus  cause  a  sizeable 
error  amounting  to  the  factor  2* 

We  shall  show  that  the  neglect  of  t-be  surface  term 
1s  Justifiable  when  momentum  Is  used  and  not  justifiable 
when  energy  is  used*  The  reason  for  this  is  that  the  rapid¬ 
ly  rising  fluid  Immediately  surrounding  the  sphere  carries 
considerable  energy  but  not  much  momentum*  This  portion  of 
the  fluid  includes  the  thin  sheath  and  the  thick  base  of  the 
surface  rise  adjacent  to  the  sphere,  all  of  which  will  bo 
grouped  together  and  called  the  splash*  We  may  thus  think 
of  the  splash  as  carrying  the  energy  which  would  normally  be 
dissipated  into  heat  in  lnelastio  impact* 

We  shall  show  that  the  percentage  of  energy  carried 
by  the  splash  is  much  larger  than  the  percentage  of  momentum 
carried  by  it,  at  least  for  the  depths  of  penetration  in 
which  we  are  interested*  This  is  plausible  on  intuitive 
grounds  since  the  mass  of  the  splash  Is  concentrated  near  the 
sides  of  the  sphere  and  is  moving  outnard  as  well  as  upward 
(see  Plate  5  )•  The  outward  or  horizontal  motion  of  the  splash 
contributes  to  the  kinetic  energy,  but  contributes  nothing  to 
the  momentum  because  of  the  rotational  symmetry  about  the  verti¬ 
cal  direction** 

At  each  instant,  the  additional  water  entering  the  splash 
is  shed  off  the  aide  of  the  sphere  at.  the  place  where  the  water 
level  intersects  the  sphere*  Let  u^  and  u^  be  the  horizontal 
and  vertical  components  of  the  splash  as  it  leaves  the  sphere, 
and  let  0  be  the  angle  indicated  In  the  diagram*  The  velocity 
with  which  the  splash  is  moving  in  the  direction  normal  to  the 
sphere  is  TJ  cos  0,  where  TJ  is  the  vortical  velocity  of  the 

*  This  remark  is  not  valid  in  the  case  of  oblique  entry* 

Because  of  the  asymmetry  there  is  then  a  sizable  horizontal 
component  of  the  momentum  of  the  splash* 
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sphere.  We  therefore  have  (see  the  diagram) 

u^  =  (u^  +  U)  cot  9 

and  so  the  magnitude  u  of  the  resultant  splash  velocity  Is 
given  by 

u  =  +  (1  +  ^)2  cot2  9  . 


Let  M  be  the  virtual  mass  of  the  water  defined  in 
Sections  3,  4*  By  (10.4),  1 =  2M  and  so  the  total  kinetic 

energy  of  the  fluid  is  ili  U2  =  MU2.  Hence,  if  in  is  the 

2  ©  9  sp 

mass  of  water  entering  the  splash,  the  fraction  of  the  energy 
carried  by  the  splash  is 


(10*6) 


Z  -  8P 


MU2 


M 


A2 

v  u 


[1  +  (1 


U  ,2 

V 


cot  9  ] 


The  total  momentum  of  the  sp3.ash  Is  directed  virtually  upward 
by  rotational  symmetry  and  is  equal  to  Dividing  by 

the  effective  momentum  MU  of  the  fluid  we  obtain 

"■sp'S,  =  # 

MU  M  '  U 


I 


(10.7) 
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for  the  fraction  of  the  momentum  carried  by  the  aplaah. 
The  ratio  of  (10*6)  to  (10*7)  is 


(10,8) 


percentage  energy  in  splash  _ 
percentage  momentum  in  splash  “ 


i  ^  [1  +  (l  +  )2cot2©]  — - — 

tv  \  sec  9-1 


for  any  •  The  following  table  exhibits  the  minimum  values 
of  the  ratio  (10*8)  for  various  angles  9: 


9 

0° 

15° 

30° 

45° 

60°  75°  90° 

minimum  possible 
value  of  ratio  (10.8) 

O0 

28.3 

6.5 

2.4 

1  .35  0 

This  table  gives  the  minimum  value  of  the  ratio  (10*8)  as  the 

fluid  is  shed  off  into  the  splash.  The  ratio  (10.8)  for  the 

total  splash  is  some  sort  of  integral  of  - 1 -  .  Thus,  in 

sec  9-1 

the  early  stages  of  the  entry  of  the  sphere  the  energy  loss 
in  the  splash  is  much  greater  than  the  momentum  loss.  For 
this  reason,  an  approximation  based  on  momentum  is  more  accurate 
than  one  based  on  energy. 

In  Sections  7,  8  of  the  present  report,  we  have  computed 
the  correction  to  the  virtual  mass  M  due  to  the  motion  of  the 
surface.  The  surface  correction  to  the  dimensionless  virtual 
mass  m  is  J' ^  Ksdb  •  This  correction  exhibits  the  behavior 
o 

indicated  above.  At  the  beginning,  near  b  =  0,  the  correction 
is  relatively  small.  The  relative  magnitude  of  the  correction 
Increases  until  near  b  =  1  it  amounts  to  about  20  percent* 
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Prom  the  mathematical  aide  we  can  aee  why  the  surface 

correction  la  smaller  on  a  momentum  basis  than  on  an  energy 

basis*  The  surface  correction,  on  a  momentum  basis,  is 

merely  dxdy.  On  an  energy  basis  the  surface  correction 

is  by  (10*5)  /°ff$  (-  )  dS  which  has  the  extra  factor  ^  . 

s  dn  3n 

The  surface  of  the  water  immediately  adjacent  to  the  sphere 

is  moving  up  rapidly,  so  that  ^  is  large  there  and  l^)dS 

/■/•.  on  '  s  °n 

is  larger  than  pJJ  $  dxdy* 

Also,  since  the  surface  is  rising,  -  ^  is  negative  (  n  is 
the  inward  pointing  normal)  and  so  dS  is  positive, 

remembering  that  is  negative  on  the  free9surface  .  Thus 


M  >  pff+ 


dxdy  <  M, 


the  inequalities  being  due  to  the  surface  terms. 
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SMOOTH  ENTRY  (Continued) 

tool  3 phene  of  diameter  2.86  cm.  dropped 
from  a  height  of  16.2  cm. 
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These  pictures  were  taken 
of  the  Morris  Dan  Group  at  California 


by  Dr.  Hurley 
Institute  cf  Technology. 
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Steel  ball,  5/8r  diace  ter,  velocity  19.4  ft/sec. 
Pressure . above  liquid  =  20  ca<  mercury. 


These  pictures  ■were  taken  by  R.  H.  Davies 
of  the  Engineering  Laboratory,  Cambridge,  England. 
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HIGH-SPEED  ENTRY  SHOWING  COMPRESSIBILITY  EFFECTS 


4 

TJ0  =  8*84  x  10  cm /sec.  A  =  0*159  cm.  Mq  =0.6  grams. 

This  photograph  was  taken  by  E.  Newton  Harvey 
of  Princeton  University. 
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